Generation of nonclassical microwave states using an artificial atom in ID open space 
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We have embedded an artificial atom, a superconducting transmon qubit, in a ID open space 
and investigated the scattering properties of an incident microwave coherent state. By studying the 
statistics of the reflected and transmitted fields, we demonstrate that the scattered states can be 
nonclassical. In particular, by measuring the second-order correlation function, </ 2 ', we show photon 
antibunching in the reflected field and superbunching in the transmitted field. We also compare 
the elastically and inelastically scattered fields using both phase-sensitive and phase-insensitive 
measurements. 

PACS numbers: 42.50.Gy,03.67.Hk,42.50.Ar,85.25.Cp 



A single atom interacting with propagating electro- 
magnetic fields in open space is a fundamental system 
of quantum optics. Strong coupling between a single 
artificial atom and resonant propagating fields has re- 
cently been achieved in a ID system experimentally 
demonstrating nearly perfect extinction of the forward 
propagating fields [2j. However, this extinction can be 
explained by classical theory: a classical point-like oscil- 
lating dipole perfectly reflects resonant incident fields [1] . 
In this Letter, we demonstrate the quantum nature of the 
transmitted and reflected fields generated from our artifi- 
cial atom in ID open space by using a resonant coherent 
state as the incident field. In particular, by measuring 
the statistics of the fields we show that the reflected field 
is antibunched [H, HJ while still maintaining first-order co- 
herence. Moreover, we observe superbunching statistics 
in the transmitted fields Q. 

To understand how our artificial atom generates anti- 
bunched and superbunched states, it is helpful to con- 
sider the incident coherent state in the photon number 
basis. For a low power incident field with less than 0.5 
average photons per lifetime of our atom, we can safely 
approximate the coherent field using only the first three 
photon eigenstates. If we consider a one-photon inci- 
dent state, the atom reflects it, leading to antibunching 
statistics in the reflected field. Together with the zero- 
photon state the reflected field still maintains first-order 
coherence. For a two-photon incident state, since the 
atom is not able to scatter more than one photon at a 
time, the pair has a much higher probability of transmis- 
sion, leading to superbunching statistics in transmission 
0, 0]. In this sense, our single artificial atom acts as 
a photon-number filter, which extracts the one-photon 
number state from a coherent state. This represents a 
novel way to generate photon correlations and nonclas- 
sical states in the microwave compared to other recent 
work 0-EI1- 

Our system consists of a superconducting transmon 




qubit |12(, strongly coupled to a ID transmission line in 
a coplanar waveguide configuration(see Fig. 1A). The 
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FIG. 1. (A) A micrograph of our artificial atom, a supercon- 
ducting transmon qubit embedded in a ID open transmission 
line. (Zoom In) Scanning-electron micrograph of the SQUID 
loop of the transmon. (B) Schematic setup for measurement 
of the second-order correlation function. This setup enables 
us to do Hanbury Brown- Twiss measurements between out- 
put ports 1 and 2. Depending on the choice of input port, 
we can measure g^ of the reflected/transmitted field. (C) 
Transmittance on resonance as a function of incident power. 
(Inset) A weak, resonant coherent state is reflected by the 
atom. 



ground state |0) and first excited state |1) have a tran- 
sition energy hu>oi. The relaxation rate of the atom is 
dominated by an intentionally strong coupling to the Z 
= 50 Q transmission line through the coupling capacitor 
C c , defined in Fig. IB. 

The electromagnetic field in the transmission line is de- 
scribed by an incoming voltage wave Vi n , a reflected wave 
Vr and a transmitted wave Vt- In Fig.lA, the transmit- 
tance is defined as T = |Vr/Vi n | 2 . For a weak coherent 
drive on resonance with the atom, we expect to see full 
reflection of the incident signal [1, [HI . This can be un- 
derstood in terms of interference between the incident 
wave and wave scattered from the atom, which destruc- 
tively interferes in transmission and constructively inter- 
feres in reflection [3, EH- In the sample measured here, we 
achieved extinction of more than 99% in transmittance, 
as shown in Fig. 1C. By measuring the transmission coef- 
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ficient as a function of probe frequency and probe power 
P, we extract w i/2tt = 5.12 GHz, T w /2ir = 41 MHz, 
T t f,/2TT = 1 MHz. The relaxation rate r 10 is dominated 
by the coupling to the line and is much greater than the 
pure dephasing in our system. We define an aver- 
age number of photons per interaction time 2tt/Tiq as 
N = 2tt P/(hoJ p T w ). N = 1 for a power of -128 dBm. 

The resonant electromagnetic field reflected from or 
transmitted through the atom (depending on the choice 
of input port, see Fig IB) is fed through two circulators 
to a commercial 90° hybrid coupler, with its other input 
terminated with 50 f2. The hybrid coupler effectively acts 
as a microwave beam splitter. Ideally, the signal coming 
into the other input of the hybrid coupler should be vac- 
uum. The two outputs of the beam splitter are sent to 
two nominally identical microwave HEMT amplifiers at 
4.2 K that have system noise temperatures of 7 K. We 
make the essential assumption in our analysis that the 
noise added by the two amplifiers is uncorrelated. After 
further amplification, the two outputs are fed into a pair 
of vector digitizers, which capture the voltage amplitude. 
We then can choose to digitally filter the voltage data 
to a desired bandwidth, BW. This setup enables us to 



perform Hanbury Brown- Twiss (HBT) [14j-|l8| measure- 
ments using linear quadrature detectors, i.e., amplifiers. 

Second-order correlation measurements provide a 
method to characterize microwave states generated in our 
system. In particular, they provide a statistical tool to 
show that the scattered light is nonclassical. For refer- 
ence, we first considered a thermal state and a coherent 
state. The thermal state [3| was generated by simply 
amplifying the noise of a 50 resistor through room 
temperature amplifiers before sending the signal down 
through the transmission line in our set-up, with the 
qubit off-resonance. The two digitizers were set to a sam- 
pling frequency of 10 8 samples/sec and then the voltages 
were digitally filtered separately. Finally, we determine 
the power-power correlations as a function of delay time 
between the two outputs. 

The second order correlation function can be expressed 

as 



(AP 1 (t)AP 2 (t + r)) 



[(P 1 (t)}-(P hN (t)}][(P2(t)}-(P2,N(t))r 



where r is the delay time between the two digitizers, 
Pi,P% are the output powers in ports 1 and 2, respec- 
tively. P\.n i P2.N are the amplifier noise in ports 1 and 
2, respectively, when the incident source is off. There- 
fore, [(Pi(t)) — (Pi,N(t))] represents the net power of the 
field from output port i, where i = 1,2. (A Pi AP 2 ) is the 
covariance of the output powers in port 1 and 2, defined 
as ((Pi — (Pi))(p2 — (P2)))- The following assumptions 
were made: (1) the amplifier noise originating from the 
two independent detection chains is uncorrelated, and 
(2) the 50 f2 terminator is in its ground state (5 GHz 
photons at ~50 mK). In Fig. 2 A, we show </ 2 ' as a func- 



tion of delay time r, for a thermal state with two differ- 
ent filter bandwidths, and also for a coherent state. For 
thermal states, g^ 2 \0) = 2 regardless of the filter band- 
width BW. The width of g^ 2 \r) for the thermal state 
is determined by the filter function. For the filter used 
g( 2 )(-r) = 1 + e - 27T BW\r\^ gQijj cur ves of the thermal 
state in Fig. 2A show this equation with no free fitting 
parameters. We had a trigger jitter of ±1 sample between 
the two digitizers. To minimize the effect of this trigger 
jitter, we oversample and then digitally filter (average) 
the data in all the g^ measurements. For a coherent 
state, we expect g^ 2 \r) = 1. This is indeed what we find 
if our atom is off-resonance from our applied coherent 
source. 

After making these initial measurements, we made sec- 
ond order correlation measurements to study the field 
transmitted through our atom. We applied an on- 
resonance microwave drive and measured g^ 2 \r) in trans- 
mission for different incident powers as shown in Fig. 
2B. Here the sampling frequency was again set to 10 8 
samples/s with BW=55 MHz. At the lowest power 
(P = —129 dBm, N = 0.8) that we can readily mea- 
sure, we see superbunching of the photons as expected 
[3, with g^ir = 0) = 2.31 ± 0.09 > 2. Superbunching 
occurs because the one-photon state of the incident field 
has been selectively filtered and reflected while the two- 
photon state is more likely transmitted (the three-photon 
components and higher are negligible). This transmitted 
state generated from our qubit is thus bunched even more 
than a thermal state. For high powers, where N ^> 1, 
we find g^ (r) = 1 as expected, because most of the co- 
herent signal then passes through the transmission line 
without interacting with the atom owing to saturation of 
the atomic response. The correlation once again looks 
like a coherent state. For all measurements shown here 
we find, g( 2 \oo) — 1, as expected. In the right inset of 
Fig. 2B, we plot g^ (0) as a function of incident power, 
and clearly see the bunching behavior decrease as the in- 
cident power increases. For comparison, we also plot the 
coherent state and thermal state. 

In Fig. 2C, we show the measured g^(r) of the re- 
flected field from our atom. At low powers, where N <C 1, 
we clearly observe antibunching of the field Each 
trace here was collected, computed and averaged over 17 
hours, corresponding to 2.4 x 10 11 measured quadrature 
field samples (2 Tbyte of data). The antibunching be- 
havior at the lowest power(P = —132 dBm, N = 0.4), 
g( 2 \0) = 0.51 ± 0.05, reveals the quantum nature of the 
field. Ideally, we would find g^(0) = as the atom can 
only absorb and emit one photon at a time. The non-zero 
(/ 2 )(0) we measured is most likely from the combination 
of trigger jitter between the two digitizers and the filter 
bandwidth (here 55 MHz). In Fig. 2D, we clearly see 
how the antibunching dip depends on bandwidth. For 
small BW, where BW <C Toi, the antibunching dip van- 
ishes entirely. With a small BW or long sampling time, 
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FIG. 2. Second-order correlation function of a thermal state, a coherent state and the states generated by the artificial atom. 
A) g^ 2 ' of a thermal state and a coherent state as a function of delay time r. B) g^ of the resonant transmitted microwaves 
as a function of delay time for four different incident powers. Inset: g (2) (0) as a function of incident power. For comparison, 
the result for a thermal state and a coherent state are also plotted. We see that the transmitted field statistics(red curve) 
approaches that of a coherent field at high incident power, as expected. For a coherent state <jr 2 '(0) = 1 is independent of 
incident power (blue). C) of a resonant reflected field as a function of delay time for two different incident powers. The 
antibunching behavior reveals the quantum nature of the field. The curves shown here had a digital filter with a 55 MHz 
bandwidth applied to each detector. D) g^ of a resonant reflected field as a function of delay time at -131dBm incident power 
for different filter bandwidths. As the bandwidth decreases, the antibunching dip vanishes. The solid curves in B), C) and D) 
are the theory curves, including trigger jitter model (see text). We extract a temperature of 50 mK from these fits, with no 
additional free fitting parameters. The peculiar feature of g^ around zero in the theory curves in B) is due to the trigger jitter 
model. The error bar indicated for each data set is the same for all the points. 



the time dynamics of antibunching cannot be resolved. 
In other words, within that time interval, the atom has 
had time to absorb and emit many photons. 

To obtain the theoretical results in Fig. 2 B, C, and 
D, we use a master equation formalism. The digital filter 
is modeled by a single-mode resonator. A master equa- 
tion describing both the transmon and the resonator is 
derived using the formalism of cascaded quantum sys- 
tems. To model the effect of the trigger jitter, the value 
of gW(r) at each point is replaced by the average value 
of 5 (2) (t), 5 ( 2 )(r-10 ns) and g (2 \T+10 ns). 

Our artifical atom selectively filters out the Fock state 
n = 1 from the input coherent state. As a result, the 
reflected and transmitted field display antibunched and 



superbunched statistics, respectively. Thus, the qubit 
acts as a passive photon-number filter, converting a co- 
herent microwave state to a nonclassical one, with high 
production rate. 

While the scattered field requires a purely quantum 
description, it can still maintain first-order coherence 
similar to a classical field, as shown below. We can de- 
fine the first-order correlation function in steady state as 
gi 1 ) = (V) 2 / (V 2 ). First-order coherence then refers to 
g^ = 1. As expected, for a thermal source, this function 
is and for a coherent state it is 1. 

The first-order coherence properties of the scattered 
resonance field strongly depend on the Rabi frequency 
fL and the relaxation rate T±a of the atom. The Rabi 
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FIG. 3. Elastic vs. inelastic scattering from the artificial 
atom. A) A microwave pump is applied at wrji- As the 
power of the a>oi pump increases, the Mollow triplet appears 
in the spectrum with peak separation equal to the Rabi fre- 
quency Q p . (inset) Dressed state picture of the energy levels. 
B) The coherently/elastically reflected power (phase-sensitive 
average, red curve) or total reflected power (phase-insensitive 
average, green, grey and blue curves) as a function of res- 
onant incident powers for different bandwidths (BW). The 
total power reflected is the sum of both the elastic and in- 
elastic fields, after we subtract off the amplifier noise power. 
Solid curves are the theory fits to experimental data. The 
black dash curve shows the input power for comparison. At 
low powers, TV << f, we observed that the reflected field 
is mostly first-order coherent: the elastically reflected power 
is the same as the total reflected power. At high powers, 
N > 1, more and more photons are inelastically scattered as 
the Mollow triplet begins to emerge. The wider the measure- 
ment bandwidth, the more of the Mollow triplet we capture. 
Note that the output power includes the 79 dB gain of the 
amplifiers. 



frequency f2 p is linearly proportional to the amplitude 
of the drive, fl p oc y/~P [2|]. For n p <C T\o, we expect 
the scattered fields to be mostly coherent to first order, 
with a spectrum well described by a delta function Q 
at wqi. For Q p > Tiq, the scattered fields contain three 
additional inelastic Lorentzians (known as the Mollow 
Triplet [2(|) centered on the frequencies u = u>oi, w = 
luqi +Cl p and u) = ojqi — Cl p , indicated in the theory plot 
Fig. 3A. 

We send a single tone at cj i and measure these scat- 
tered (reflected) fields from only one of the output ports, 
as shown in Fig. 3B. Note that we see the same behav- 
ior from both outputs. We use a phase-sensitive average 



(V) to capture the elastic (coherent) component of the 
scattered field. For the total scattered field, the sum of 
the elastically and inelastically scattered fields, we use a 
phase- insensitive average (V 2 }. The amount of the in- 
elastic field that we capture depends on the bandwidth 
(BW) of our measurement, as indicated in Fig. 3A. The 
solid curves are the theory fits, using the model in Fig. 3A 
(integrating the Mollow triplet), with no free fitting pa- 
rameters (using the same parameters extracted before). 
As expected, at low incident power, the total scattered 
field is roughly equal to the elastically scattered (coher- 
ent) field. This indicates that the field is first-order co- 
herent with ~ 1. We note that this is also the regime 
where antibunching is observed. At high incident fields, 
where il p > Tiq, the main contribution to the total field 
is from inelastic scattering. 

In conclusion, we investigated the scattering properties 
of a single artificial atom in a ID open space by measuring 
the first, second and fourth moments of the voltage field. 
We verified the quantum nature of the scattered field 
using the second-order correlation function, while also 
showing that the field maintained first-order coherence. 
In fact, the whole process leads to a redistribution of the 
photon number state @. 

We plan to investigate applications of this phe- 
nomenon, such as generating single-photon states on de- 
mand. This system may offer advantages over placing an 
artificial atom in a cavity [21 - 23 1 , for instance. In our 
system, regardless of the strength of the input pulse only 
a single photon state is reflected. This eliminates the 
need for precise Rabi pulses to excite the atom as in a 
cavity-based photon source. Furthermore, our generated 
single photons can maintain the same envelope as the 
incident coherent state with a wide bandwidth limited 
by the atom relaxation rate. For a cavity-based photon 
source, the bandwidth is limited by the cavity width and 
subject to the problem of stochastic release by the cavity. 
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